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(54) Title: DNA MARKERS FOR MEAT TENDERNESS 

(57) Abstract: A method Tor assessing the tenderness of meat from an animal, comprising the step of testing ihe animal for the 
presence or absence of a genetic marker selected from the group consisting of: ( I) an allele of the gene encoding caJpastatin (CAST) 
associated with peak- force variation or genetic variation located other than in the CAST gene which shows allelic association with 
the CAST allele; and (2) an allele of the gene encoding lysyi oxidase (LOX) associated with intron compression of the semitendinous 
muscle or genetic variation located other than in the LOX gene which shows allelic association with the LOX allele. 
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DNA MARKERS FOR MEAT TENDERNESS 
Technical Field 

The present invention is concerned with genetic 
markers for meat tenderness in animals, and with methods 
5 and oligonucleotide probes for assessing meat tenderness 
in said animals, and a kit for this purpose. The 
invention is useful for the selection of animals which 
show desirable traits in meat tenderness either for 
breeding or to select animals destined to be slaughtered 
10 for food. 

Background Art 

Meat tenderness is an important issue for 
consumers, and one which can influence demand sufficiently 

15 for an especially tender meat to command a premium price 
in the marketplace. The physiological change in muscle 
structure during the postmortem period is complex but 
clearly seems to be at least one factor in meat 
tenderness. The calpain/calpastatin system is an 

20 endogenous, calcium-dependent proteinase system, theorised 
to initiate in vivo muscle protein degradation. 
Calpastatin appears to inhibit calpain activity and 
therefore may be assumed to have a role in meat tenderness 
through the regulation of postmortem proteolysis. In 

25 particular, calpain is response for the breakdown of 
myofibril protein, which is closely related to meat 
tenderness. 

It might therefore be suspected that calpastatin 
activity will be related to meat tenderness. Indeed, an 

3 0 increase in postmortem calpastatin activity has been 
correlated to reduced meat tenderness. Nevertheless, 
despite such observations, no clear link between the CAST 
gene, which encodes calpastatin, and meat tenderness has 
been established. 

3 5 For example, Lonergan et al. (1995) undertook a 

restriction fragment length polymorphism analysis at CAST 
and failed to find an association with either calpastatin 
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activity or tenderness in cross bred offspring of sires 
from eight breeds. Chung et al. (1999) measured 
calpastatin activity, Warner-Brat zler Shear Force and 
myofibril fragmentation index in forty-seven purebred 
5 Angus bulls. However, they concluded that U PCR single- 
strand conformation polymorphism analysis of the 
calpastatin gene was not useful for prediction of 
calpastatin activity, myofibril fragmentation index or 
meat tenderness". 

10 It is long known that one of the actions of lysyl 

oxidase (LOX) is to initiate crosslink formation at an 
early stage in collagen f ibrillogenesis (e.g., Cronlund et 
al., 1985). The action of lysyl oxidase is intensively 
studied with hundreds of publications on a variety of 

15 aspects of its importance in cancer (Giampuzzi et al., 
2001), the vasculature (Nellaiappan et al.) and other 
tissue and organ systems. 

Variation at the gene itself has not been 
associated with differences in beef tenderness although 

20 LOX has always been seen as a strong candidate on 
biochemical grounds for a gene contributing to the 
collagen component of tenderness. Analysis of genetic 
linkage has implicated the genomic region that includes 
LOX in linkage analysis of family variation in adhesion 

25 and instron compression of the semitendinosis muscle 

(STADH and STIC; Drinkwater et al., 1999). However, LOX 
itself has not been associated with these measures of 
tenderness through the study of population associations. 

Meat tenderness is a complicated trait because 

3 0 there are many sources of variation that affect postmortem 
meat tenderisation. Some of these are non-genetic effects 
such as the age of the beast, the nature of its feed, 
degree of stress prior to slaughter, carcass chilling, 
postmortem ageing time and cooking and testing methods. 

35 It has been suggested (e.g. Koohmaraie (1994)) that 

approximately 30% of the variation in tenderness in meat 
can be explained by additive gene effects within a single 
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breed, and that approximately 10% of the variation is 
explained by environmental and non-additive gene effects. 
In the Lonergan study the cattle were slaughtered at just 
over 1 year of age (430 days), the sample contained only 
5 83 animals of random peak-force values, and the sample 
consisted entirely of crosses between various taurine 
breeds. Likewise, in the Chung study purebred Angus bulls 
only 280 days of age were used. In addition, in neither 
study were the animals selected for extreme peak-force 
10 values, and it therefore seems that environmental and. non- 
fixed genetic effects may have contributed to the failure 
to identify any genetic linkage between the CAST gene and 
meat tenderness. 

15 Summary of the Invention [Revise this once claims are 
settled] 

Through using a protocol designed to reduce the 
influence of fixed genetic and environmental effects, the 
present inventor was unexpectedly able to show allelic 

20 association between the CAST and LOX genes and meat 

tenderness. In general terms, therefore, the present 
invention is concerned with genetic markers for meat 
tenderness in animals killed for meat which are 
polymorphisms of the CAST and LOX genes or polymorphisms 

25 which show allelic association therewith. 

Accordingly, in a first aspect of the present 
invention there is provided a method for assessing the 
tenderness of meat from an animal, comprising the step of 
testing the animal for the presence or absence of a 

3 0 genetic marker selected from the group consisting of: 

(1) an allele of the gene encoding calpastatin 
(CAST) associated with peak-force variation or 
genetic variation located other than in the CAST 
gene which shows allelic association with the 

3 5 CAST allele; and 

(2) an allele of the gene encoding lysyl oxidase 
(LOX) associated with variation in instron 
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compression of the semitendinosis muscle or 
genetic variation located other than in the LOX 
gene which shows allelic association with the 
LOX allege. 

5 Preferably, the allele tested for is located in 

the 3' UTR of CAST, and is typically CAST 3 D/E allele 1, 
having the following partial DMA sequence: 

catttggaaaacgatgcctcacgtgttcttcagtgttctgatttctcat 
gacccctttcctcttGgacttgtgggactgtgtttgatgtttccctgggttgttgttt 
10 ataagtcagtcataaAatactgtgcattgggcacatgtctcctcttgagctgctaatc 
gtaga (SEQ ID NO:l) , 

CAST3 D/E allele 2, having the following partial 
DNA sequence: 

catttggaaaacgatgcctcacgtgttcttcagtgttctgatttctcat 
15 gacccctttcctcttAgacttgtgggactgtgtttgatgtttccctgggttgttgttt 
ataagtcagtcataaAatactgtgcattgggcacatgtctcctcttgagctgctaatc 
gtaga (SEQ ID NO: 2) 

or CAST3 D/E allele 3, having the following 
partial DNA sequence: 
20 cat ttggaaaacgatgcctcacgtgttcttcagtgttctgatttc teat 

gacccctttcctcttAgacttgtgggactgtgtttgatgtttccctgggttgttgttt 
ataagtcagtcataaTatactgtgcattgggcacatgtctcctcttgagctgctaatc 
gtaga ( SEQ ID NO : 3 ) * 

Reduced meat toughness is selected for by 
25 rejecting animals with the "11" and u 12" genotypes and 

accepting animals with the "22" or "23" genotypes. In the 
sequences given above, the allelic difference is 
highlighted with a capital letter. These three alleles in 
the D/E DNA fragment are due to two SNP (single nucleotide 
30 polymorphisms). The first SNP is at base 2655 of Genbank 
sequence L14450, which is the same as base 2959 of Genbank 
sequence AF159246; it is a G to A change so that allele 1 
has G and alleles 2 and 3 have A. The second SNP is an A 
to T change 58 base pairs 3' to the first SNP. Since only 
3 5 three alleles have been noted for this region, with 2 
' SNPs, it implies that there is complete linkage 

disequilibrium between allele 2 and allele 3, and allele 
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3 may have evolved from allele 2. This is expected since 
they are 58 base pairs apart. For predictive purposes, a 
test of the second SNP which gives a result of allele 3 is 
equivalent to a test of the first SNP giving a result of 
5 allele 2. This is consistent with the peak force values 
of animals that are '23' heterozygotes, all of whom have 
low peak force values- While not wishing to be bound by 
theory, it is believed that these polymorphisms are linked 
to a mutation in or near the calpastatin gene (perhaps in 

10 the promoter or an intron) which results in reduced 
calpastatin expression or activity. 

A further polymorphism has been identified in the 
5' UTR of the CAST gene and other polymorphisms which 
exhibit allelic association with the polymorphism of the 

15 3' UTR, and therefore also act as genetic markers for the 
tenderness characteristics described above, may also be 
present at least within the genomic DNA embracing the 
coding region of the CAST gene and the 5' UTR and 3' UTR 
regions of that gene. In addition, where there has been a 

20 recent reduction in population size for a species, 

particular haplotypes of individuals will be relatively 
over-represented. If insufficient time has elapsed to 
cause allelic association to decay, there will be linkage 
disequilibrium even for alleles which are far apart. 

25 Livestock species such as cattle have been domesticated 
from a relatively small pool of wild ancestors in recent 
times, and therefore in these species allelic association 
is found between alleles that may be remote physically. 
Thus, it may be expected that regions of genetic variation 

30 that are outside the CAST gene will also show allelic 
association with the polymorphisms in the CAST gene 
described above, and therefore will be suitable genetic 
markers for the characteristic of peak-force variation. 
Hence, these polymorphisms may also be used to assess meat 

35 tenderness. 

In particular the CASTS microsatellite 
polymorphism (Nonneman et al, 1999) has been found to be 
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useful as a genetic marker for meat tenderness. Of the 
more common alleles, alleles 7 and 9 have been found to be 
associated with low peak-force and allele 3 to be 
associated with hi.gh peak-force. 
5 Therefore, the invention encompasses, in 

preferred embodiments, the further step of testing for the 
presence or absence of one or more additional genetic 
markers such as alleles of the gene encoding calpastatin 
associated with peak-force variation, in particular, with 

10 testing for the presence or absence of CASTS allele 7 or 9 
and/or the presence or absence of CASTS allele 3. The 
most favorable results when the presence of CAST D/E 
allele 2 has been established is to have CASTS allele 7 or 
allele 9 present also, whereas the benefits of the 

15 presence of CAST D/E allele 2 are to some degree 

counteracted if the animal also possesses CAST5 allele 3. 

The LOX polymorphism has also been shown to be a 
genetic marker for meat tenderness, and allele 1 or allele 
2 may be tested for. Just as for the CAST gene, allelic 

20 association may be exhibited to alleles located outside 
the LOX gene. 

According to a second aspect of the present 
invention, there is provided a genetic marker for meat 
tenderness in an animal which is a polymorphic form of the 

2 5 CAST gene, being the CAST3 D/E polymorphism or the LOX 

polymorphism. 

According to a third aspect of the present 
invention there is provided an isolated DNA molecule 
comprising the nucleotide sequence set forth in SEQ ID 
30 NO:l, SEQ ID NO: 2 or SEQ ID NO: 3. 

According to a fourth aspect of the present 
invention there is provided an isolated DNA molecule 
consisting of the nucleotide sequence set forth in SEQ ID 
NO:l, SEQ ID NO:2 or SEQ ID NO: 3. 

3 5 According to a fifth aspect of the present invention 

there is provided a method for selecting an animal likely 
to yield meat of improved tenderness, comprising the steps 
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of: 

(1) testing the animal for the presence of an allele 
of the gene encoding calpastatin (CAST) 
associated with low peak-force or genetic 

5 variation located other than in the CAST gene 

which shows allelic association with the CAST 
allele and/or for the presence of an allele of 
the LOX gene associated with the low instron 
compression of the semitendinosis muscle or 
10 genetic variation located other than in the LOX 

gene which shows allelic association with the 
LOX allele; and 

(2) selecting animals which have the CAST and/or LOX 
allele and/or genetic variation in allelic 

15 association therewith. 

Advantageous ly, in order to assess the tenderness 
of meat from an animal and/or select an animal likely to 
yield meat of improved tenderness testing may comprise the 
steps of: 

20 (1) obtaining a biological sample from the animal; 

(2) extracting DNA from the sample; 

(3) amplifying DNA from the CAST or LOX gene and/or 
from regions of genetic variation which show 
allelic association to polymorphisms of the 

25 relevant one of the CAST or LOX gene; and 

(4) identifying the allele present in the amplified 
DNA. 

Typically the allele identified in step (4) is 
one of CAST 3 D/E allele 1, CAST 3 D/E allele 2 and CAST 3 
3 0 D/E allele 3 described above and/or the CAST 5 alleles 
described above. 

Preferably the biological sample is blood, but 
other biological samples from which DNA can be amplified 
may be used. For example, hair root samples, cheek 
3 5 scrapings, skin samples and the like may be used. 

Typically amplification is performed using the 
polymerase chain reaction (PCR) , but other DNA 
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amplification methods such as the ligase chain reaction 
are well known in the art, and may alternatively be used. 

Preferably the alleles are identified by 
polyacrylamide gel electrophoresis techniques such as 
5 SSCP, or by other techniques well known to the person 
skilled in the art such as RFLP analysis. 

In a sixth aspect the invention provides an 
oligonucleotide probe for amplification of a genetic 
marker associated with peak-force variation, said genetic 
10 marker being either an allele of the gene encoding 

calpastatin (CAST) or genetic variation located other than 
in the CAST gene which shows allelic association with said 
allele. 

Typically the probe is selected from the group 
15 consisting of: 

castd 5' cat ttg gaa aac gat gcc tea c 3' 
caste 5' tct acg att age age tea aga gga g 3' 
CAST5U1 5 ' -GTAAAGCCGCACAAAACACACCCAGG- 3 ' 
CAST5D1 5 ' -GTTTCTGGACCCTCTGGATGAGGAAGCGG- 3 ' . 
20 In view of the designation of the primers as 

CASTD and CASTE, the amplified fragment of the CAST gene 
is referred to from time to time as the CAST D/E fragment 
and the polymorphism as the CAST D/E polymorphism. 

According to a seventh aspect of the present 
25 invention there is provided an oligonucleotide probe for 
amplification of a genetic marker associated with 
variation in instron compression of the semitendinosis 
muscle, the genetic marker being either an allele of the 
gene encoding lysyl oxidase (LOX) or genetic variation 
3 0 located other than in the LOX gene which shows allelic 
association with said allele. 

Typically the oligonucleotide probe is an 
oligonucleotide probe selected from the group consisting 
of: 

3 5 LOX K5 : 5' tat cac tga tgt caa acc tg 3' 

LOX K6: 5' act cag gca cca aat age tg 3' 
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According to an eighth aspect of the present 
invention there is provided a kit for use in assessing the 
tenderness of meat from an animal and/or selecting an 
animal likely to yield meat of improved tenderness, 
5 comprising oligonucleotide probes for amplification of at 
least one genetic marker for meat tenderness, said genetic 
marker being either an allele of the gene encoding 
calpastatin (CAST) or genetic variation located other than 
in the CAST gene which shows allelic association with said 

10 allele, or an allele of the LOX gene associated with low 
instron compression of the semitendinosis muscle or 
genetic variation located other than in the LOX gene which 
shows allelic association with the LOX allele, and means 
for amplifying DNA. 

15 The primers used to amplify the DNA are the CASTD 

and CASTE primers and/or the CAST5U1 and CAST5D1 primers 
for amplifying the CASTS polymorphism. However, other 
primers able to amplify polymorphisms associated with a 
reduction in toughness in meat are envisaged, whether 

20 these be primers which amplify a polymorphism other than 

the CAST 3 D/E polymorphism or CASTS polymorphism, or other 
primers able to amplify the CAST 3 D/E fragment of CAST 5 
polymorphism. 

The methods of the invention may be used both for 

25 the selection of breeding animals and for the selection of 
unpedigreed animals for entry into feed lots. In the 
latter case, the methods of the invention allow for 
animals with unsuitable pedigrees to be excluded from feed 
lots on the basis that highly tender meat is unlikely to 

30 be attained with these animals even after a. long feed lot 
holdings. Alternatively, such measurements may allow for 
determination of the optimum time to reach maximum meat 
tenderness. The invention is therefore also concerned 
with animals when selected by the method of the invention, 

35 their progeny and the use of both selected animals and 
their progeny for breeding, as well as meat from these 
animals. 
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The methods of the invention are applicable to 
animals including but not limited to cattle and other 
bovids, including water buffalo and bison, to other 
ungulates, including sheep, goats and deer, and pigs and 
5 chickens . 

Throughout this specification, the words 
"comprise", "comprises" and "comprising" are used in a 
non-exclusive sense, except where the context requires 
otherwise . 

10 It will be clearly understood that, although a 

number of prior art publications are referred to herein, 
this reference does not constitute an admission that any 
of these documents forms part of the common general 
knowledge in the art, in Australia or in any other 

15 country. 

Brief Description of the Drawings 

Preferred embodiments of the present invention 
will now be described, by way of example only, with 

20 reference to the accompanying drawings, in which: 
Fig 1 is a photograph of a single strand 
conformational polymorphism (SSCP) gel which shows 
genotypes for the CAST3 D/E polymorphism, from left to 
right, 11, 22, blank, 11, 12, 12, 12, 12, 11, 22. 

25 Fig 2 a & b show the distribution of Warner- 

Bratzler peak-force measurements in the two samples of 169 
and 77 animals respectively. Note that extremes were used 
so there is no middle to the distribution. It does not 
imply that the distribution is bi-modal. Note different 

3 0 scales in the figures. 

Fig 3 a & b are a plot of the raw Wamer-Bratzler 
peak- force measurements against the CAST genotypes. Note 
the gap in the middle due to the use of extremes of the 
distribution. Note the similarity between the 

3 5 distributions in the two samples. 

Fig 4 is a boxplot of the residual Wamer- 
Bratzler peak-force measurements (XI) for each genotype 
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for the first sample. The median quarter and three- 
quarter percentiles, whiskers and outliers are shown. 

Fig 5 shows the distribution of DNA fragment 
sizes for the CASTJ5 microsatellite . Horizontal axis is 
5 the frequency of each allele and the vertical axis is the 
DNA fragment size. Alleles are labelled in increasing DNA 
fragment size so allele ml in this distribution is less 
than 132 bp, m2 < 13 6 bp, m3 < 13 8 bp, m4 < 142 bp, m5 < 
144 bp, m6 < 146 bp, m7 < 148 bp, m8 < 151 bp, m9 < 153 
. 10 bp, mlO < 155 bp, mil < 157 bp, ml2 < 159 bp, m!3 < 161 
bp. DNA fragments were not found in some of the 2 bp • 
bins, e.g., in the less than 134 bp bin, and these are 
either extremely rare or non-existent. 

Fig 6 is a box plot of raw LD peak- force values 
15 along the horizontal axis versus CAST5 microsatelite 
allele identity along the vertical axis. The boxes 
contain the median value, represented by the dot, a box 
representing the 25 and 75 percentile and whiskers 
indicating the expected range for the distributions, with 

2 0 outliers indicated by open circles. Care must be taken in 

interpreting this figure since there are some alleles that 
are rare, such as ml, m2 and ml3 (see Figure 5 for the 
full distribution, so interpretations made on those 
alleles are not particularly informative. Note 

25 particularly that this is a distribution of extreme 

values, so the median value will swing from a low to a 
high value if half the samples are high values. 

Figure 7 is a photograph of a single strand 
conformational polymorphism gel showing the genotypes of 

30 LOX from left to right, 11, 11, 12, 12, 22, 11. 

Figure 8 shows the distribution of instron 
compression measurements for the two samples of 166 and 87 
animals combined. Note that extremes were used so there 
is no middle to the distribution. It does not imply that 

3 5 the distribution is bi -modal. 

Figure 9 shows the distribution of adhesion 
measurements for the two samples of 166 and 87 animals 



i 
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combined- Although the sample was selected for extremes 
of instron compression, this has not been translated into 
a series of extreme adhesion measurements. 

Figure 10 is a plot of STIC versus STADH for the 
5 combined sample. Note the non-uniformity caused by. the 
selected STIC values. 

Modes for Performing the Invention 

Example 1 - Identification of CAST3 D/E polymorphism 
10 Cattle were chosen from the DNA Bank of the 

Cattle and Beef Cooperative Research Centre located in 
Brisbane, Australia to have as diverse a genetic and 
phenotypic background as possible. Information stored in 
the CRC Database was used to select animals. Animals of 
15 extremes of peak-force were selected, although animals 

with peak-force measures above 12 were excluded since they 
might have confounded peak-force measurements. In 
essence, the procedure was to select cattle in each 
contemporary group which were of phenotypic extreme 
2 0 measures, to ensure that no sire was represented by a 
cluster of offspring, that all markets and finishing 
regimes were included in each extreme, so that extremes 
were not biased by being representative of a particular 
market or finishing regime. A total of 169 samples were 

2 5 obtained (Table 1) for the first sample. A second sample 

of 77 animals (Table 6} were analysed as a check to 
determine whether the same allelic association could be 
observed in another sample. 

These DNA samples were genotyped for the CAST 
30 (calpastatin) D/E DNA fragment using the primers 
castd 5' cat ttg gaa aac gat gcc tea c 3' 

caste 5' tct acg att age age tea aga gga g 3'. 

The conditions of the polymerase chain reaction 
(PCR) are an annealing temperature of 60 Celsius, 2.5 mM 

3 5 Magnesium chloride, and reagent mixes obtained from 

Biotech International. The DNA fragments were labelled 
via the incorporation of 32 P dCTP into the fragments during 
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the PCR, and the fragments were visualised by 
autoradiography using X-Ray film exposed overnight at room 
temperature. Alleles were scored in numerical order where 
the fastest migrating allele is number 1. 
5 The genotypes were analysed using generalised 

linear models (GLM) following the equation peak-force = 1 
+ genotypes nested within fixed effects + error 
implemented via the S-PLUS software. Fixed effects that 
were considered were breed/ finish (Australia, Korea, 

10 Japan), contemporary group (cohort), region (pasture v 

grain, north v south) and the covariate of final weight. 
The genotypes were nested within region and breed since 
pure-bred offspring of taurine sires were not pastured in 
the north. The size of the effects associated with 

15. genotype was estimated by the comparison of variances (eg, 
Andersson-Eklund and Rendel, 1993) . To estimate the size 
of effect associated with genotypic substitution, the same 
model was fitted without the calpastatin genotypes. 
Residuals were extracted and compared to the calpastatin 

2 0 genotypes. These were analysed using an analysis of 

variance to obtain adjusted means for each genotype Plots 
of raw and residual peak-force values against calpastatin 
genotypes were constructed. 

2 5 Example 2 - Analysis of CAST 3 D/E polymorphism 

There are two common alleles (Figure 1) and at 
least one rare allele for the CAST D/E polymorphism and 
both the common alleles are found in all the breeds, 
although there are clear differences in genotype frequency 

3 0 within the breeds. Zebu breeds have a greater frequency of 

the % 11' genotype (Tables 2 and 7) than taurine breeds in 
this sample. 

The raw values (Figures 2a & 2b) were then 
plotted against the CAST genotypes (Figures 3a & 3b) and 
3 5 these associations are sufficiently strong to show visual 
associations between peak-force and genotype. The most 
important genetic effect considered in the literature for 
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CAST, breed or taurine versus zebu, has been carefully 
matched so that there are animals of high and of low peak- 
force from each breed in the sample, and breed is not 
expected to be an .explanatory variable here. 
5 The analysis (Table 3) of the CAST genotypes 

showis strong, confirmatory evidence of effects of the CAST 
gene or sequences near the CAST gene on peak-force. The 
analysis shows no effect of breed, but since the sample 
consists of individuals of high and low peak-force for 

10 each breed, this was not unexpected. The size of effect 
associated with this polymorphism is approximately 7.9 
percent of the phenotypic variance estimated as a main 
effect, and the deviance associated with CAST genotype 
nested within breed within region is 121.4 (17 df, P = 

15 0.001894). An un-nested interaction term between breed 
and CAST genotype was calculated for this sample, but is 
was not statistically significant. The GLM of the CAST 
genotypes (Table 4) against the residual peak force 
measurements show a statistically significant level of 

2 0 association similar to that of the CAST genotypes 

considered as a main effect (Table 3) rather than when 
they are nested within region and breed. 

A boxplot of genotypes versus the residual peak 
force measurements (Figure 4) shows clear differences in 

2 5 distributions and the difference between medians of the 

and the x 22' genotypes are approximately 1.2 kg of 
adjusted peak force. The adjusted means from the analysis 
of variance (Table 5) gives a difference of 1.34 kg of 
peak force between the homozygote genotypes. The overall 

3 0 standard deviation for the residuals is 1.61. 

The GLM of the confirmatory sample of 77 animals 
showed a statistically significant association between 
CAST genotype and peak force, with the *1' allele 
associated with higher peak forces. When the full model 
35 was calculated, none of the factors were statistically 

significant, possibly as a result of the relatively small 
sample size. Terms in the model were dropped one by one 
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using the reduction in AIC as the criterion. All terms 
except the calpastatin genotypes were dropped (Table 8) in 
this automatic procedure, and these show a deviance of 
17.9 (2 df, P < 0^05) explaining 9.5 percent of the 
5 phenotypic variance. This is similar to the 8.9 percent 
found when the CAST genotypes were compared without other 
factors to peak force in the first sample. 

Discussion 

10 The results presented here indicate that genetic 

variation at the CAST gene is important in explaining 
variation for Warner-Bratzler peak-force measurements 
between individuals irrespective of the breed of origin. 
The sample was chosen to control for the effects of breed 

15 and to spread the sample as widely as possible over 

different sire lines, contemporary groups, feeding and 
finishing regimes; care was taken to ensure that, as much 
as possible, individuals in either extreme were from each 
breed, contemporary group, feeding and finishing system. 

20 In this way, systematic effects of these factors on peak- 
force were controlled so that the effect of the alleles 
would not be due to inadvertently being carried along by 
other factors affecting peak-force values. Indeed, there 
are statistically significant deviations in peak-force due 

25 to allelic substitution at this locus even when there is 
no accounting for the other fixed effects. Inspection of 
the raw data show frequency differences within breeds for 
the different genotypes so that the *1 # allele is rarer in 
the extreme with lower peak-force values. 

3 0 This association between the % 1 # allele and 

higher peak force measurements is confirmed in a second 
smaller sample of extreme animals. These animals are less 
extreme than those in then first sample, they are the 
left-over extremes, and they clearly show not only that 

3 5 the calpastatin genotypes are important but that in such a 
small sample, other factors known to be important are not 
found to be statistically significant. In a well matched 
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sample such as this it is not of concern, since we attempt 
to remove the effects of the other factors as much as 
possible through the choice of samples to analyse. 

The size, of the homozygote substitution is 
5 approximately 1.34 kg of peak force for the LD, equivalent 
to 0.83 of standard deviation. This value is likely to be 
overestimated since the extremes of the distribution were 
used, and a proper estimate will require animals chosen at 
random from the full distribution of peak force. 

10 Nevertheless,- this is a useful amount of genetic variance 
associated with a single marker and it is expected that 
this marker would be useful in direct DNA marker tests for 
breeding and feedlot streaming. 

The analysis shows no evidence of a breed by 

15 genotype interaction on peak- force, which means that there 
is no evidence that the allele association is different or 
absent in some breeds. This is interpreted to mean that 
there is -no heterogeneity in the breeds for the 
association between calpastatin and peak-force. 

20 A positive test for allelic association generally 

means that the causative mutation is close to the DNA 
markers. Associations in other studies have indicated 
that allelic association decays at an extremely rapid rate 
so that DNA markers even relatively close to a 

2 5 quantitative trait locus will find no evidence of 

association (e.g., Coleman et al., 1995; Barendse, 1997). 
This indicates that the causative mutation or mutations 
are extremely close to the CAST D/E polymorphism. 

3 0 Example 3 - Identification of CAST 5 microstatellite 

' polymorphism 

To determine whether other polymorphisms in the 
CAST gene are associated with tenderness, both of the 
cattle samples (Tables 1 & 6) were genotyped with the 
3 5 CASTS microsatellite polymorphism (Nonneman et al, 1999) 
which was developed from DNA sequence reported earlier 
(Cong et al . , 1998) . 
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The primer sequences to amplify this polymorphism 

are 

CAST5U1 : 5 ' - GT AAAGCC GC AC AAAAC AC AC C C AGG - 3 ' and 
CAST5D1 :., 5' - GTTTC TGGAC C C TC TGGATGAGGAAGC GG - 3 ' 
5 with amplification fragments in the range 130 - 159 bp, 

sizes determined on an ABI 37 3 DNA sequencer. Alleles and 
genotypes were assigned based on these size fragments 
leading to 13 alleles and the distribution of allele sizes 
is shown in Figure 5, 

10 Two different sets of analyses were performed. 

In the first, the genotypes at the CAST3 D/E polymorphism 
were compared to the CASTS microsatellite to determine 
whether there' were significant associations between the 
genotypes, as a consequence of haplotypes existing along 

15 the DNA sequence. If CASTS and CAST 3 show significant 
haplotypes , since they are on either side of the CAST 
coding sequence, then all polymorphisms for the CAST 
coding sequence are expected to be in linkage 
disequilibrium with either or both of these DNA markers. 

2 0 In the second, the CASTS microsatellite alleles were 

compared to the LD peak- force measurements to determine 
whether there was any association between CASTS and 
tenderness. 

2 5 Haplotypes between CASTS and CAST 3 

Since genotypes of parents of these animals were 
not available haplotypes where determined by analysing 
animals in which one or both of CASTS and CAST 3 had 
homozygous genotypes . The frequency of these haplotypes 

30 were tabulated (Table 9). These frequencies were tested 
for heterogenity using a generalised linear model and 
found to be highly heterogenous (Table 10). This means 
that each allele at CAST 3 D/E is preferentially associated 
with specific alleles at the CASTS microsatellite. 

35 

Association between tenderness and CASTS 

Since CAST5 has 13 alleles and hence there are 91 
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possible genotypes, not all of these genotypes will be 
seen in a sample of 240 samples, as in this study, so the 
association was estimated on the alleles. As for the 
CAST3 D/E DNA marker, the polymorphism was compared to 
5 the raw LD peak-force values (Table 11 a), was examined 
for differences in interactions between breeds (Table 11 
b) , and was compared to the LD peak-force values after 
market, cohort, breed and finish effects were accounted 
for (Table 11 c) . In the last of these analyses, CASTS 

10 alleles are nested within finish and breed, as in the 
analysis of CAST3 . 

These analyses show that there is no interaction 
between CASTS allele frequency and breed on LD peak-force, 
that the association between CAST5 and the raw LD peak- 

15 force values is statistically significant at the threshold 
P < 0.01, but .when the CASTS alleles are nested within 
breed and finish, the association has a deviation which is 
0,1 > P > 0.05. The lack of interaction between CASTS and 
breed in explaining LD peak- force means that any 

2 0 differences in gene frequencies between breeds are not 

responsible for the association between CASTS and LD peak- 
force. The association between CASTS and LD peak-force in 
sections a and b of Table 11 indicate that there is some 
evidence for CASTS associated with LD peak-force, but a 
25 bias might still exist, which is why the factors such as 
market, cohort, breed and finish are also corrected for. 
Once those factors are corrected, there is a lack of 
strength in the association. In the CAST3 D/E analysis, 
correcting the additional factors improved the evidence 

3 0 for the association, and since the same samples are used, 

we know in which direction the deviations should go. Thus 
the lack of strength probably means that the large number 
of alleles nested within breed and finish, has failed to 
find an association due to the creation of a large number 
3 5 of categories. Larger numbers of alleles are expected to 
reduce the strength of associations purely due to the 
number of categories (cf Terwilliger, 1995). 
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The CASTS polymorphism can be used in conjunction 
with the CAST3 D/E polymorphism to predict LD peak-force. 
For CAST3 D/E the ell genotype is associated with higher 
peak force values , v the cl2 genotype is intermediate and 
5 the c22 genotype has the lower peak force values. 

Secondly, there is linkage disequillibrium between CAST 3 
D/E and CAST 5 . By examining the table of haplotypes, 
looking at the common microsatellite alleles, CAST 3 D/E al 
(allele 1) is most often associated with CASTS m3 (allele 

10 3) with low abundances for m7 and m9. On the contrary, 
CAST3 D/E a2 (allele 2) is most often associated with 
CASTS m9, with a similar large association to m3 and a 
lesser but still significant association with m7 . 
Inspection of Figure 6, a plot of raw LD peak-force values 

15 for each CASTS microsatellite allele, indicates that CASTS 
m7 and m9 have lower peak force values while CAST 5 m3 has 
higher peak force values. Since most of the m3 alleles 
are actually associated with CAST 3 D/E a2 and not al (108 
versus 14), this higher value is not likely to be the 

2 0 effect of CAST D/E al. Rather it provides a tool to 

refine the assignment of animals to groups, since animals 
selected for having CAST 3 D/E a2, so that they might have 
lower peak force values, might still have higher peak 
force values if they possessed CASTS m3 . They are expected 
25 to have a greater likelihood of having lower peak force 
values if they possessed both CAST 3 D/E a2 as well as 
CASTS m7 or m9 . 

Example 4 

3 0 This example shows the testing of a DNA marker in the 

LOX gene for population associations to STIC and STADH. 
Repeated statistically significant positive associations 
were found between genotypes and both STIC and STADH. 
These indicate that, unusually, the heterozygote may be 
35 one of the extreme genotypes suggesting some 

overdominance. These associations are found in a study of 
6 breeds of cattle with a structure to detect linkage 
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disequilibrium and would indicate that the gene LOX either 
contained or was located near to the genetic factor 
associated with connective tissue strength. 

Materials and Methods 

5 Cattle were chosen from the CRC DNA Bank to have 

as diverse a genetic and phenotypic background as 
possible. Two groups of animals were chosen, the first 
and larger set to test for associations and the second 
smaller set to confirm the polarity of the associations 

10 (cf. Barendse 1997; Barendse et al . , 2000). Information 
stored in the CRC Database was used to select animals. 
Animals of extremes of instron compression in the 
semitendinosis muscle were selected. Adhesion measures 
for these animals were also extracted from the database . 

15 In essence, the procedure was to select cattle in each 
cohort which were of phenotypic extreme measures, to 
ensure that no sire was represented by a cluster of 
offspring, that all markets and finishing regimes were 
included in each extreme, so that extremes were not biased 

2 0 by being representative of a particular market or 

finishing regime. A total of 253 individuals were 
selected comprising a first sample of 166 animals and a 
second sample of 87 animals (Table 11) . 

The DNA was genotyped for the LOX (Lysyl Oxidase) 

25 DNA fragment using the primers LOX K5 : 5' tat cac tga tgt 
caa acc tg 3 ' and LOX K6 : 5 ' act cag gca cca aat age tg 
3'. The conditions of the polymerase chain reaction (PCR) 
are an annealing temperature of 60 Celsius, 2.5 mM 
Magnesium chloride, and reagent mixes obtained from 

30 Biotech International. The DNA fragments were labelled 

via the incorporation of 32 P dCTP into the fragments during 
the PCR. The fragments were digested with Hinfl overnight 
at 37 Celsius before separation on gels. The fragments 
were visualised via autoradiography to X-Ray film 

3 5 overnight at room temperature. Alleles were scored in 
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numerical order where the fastest migrating allele is 
number 1 . 

The genotypes were analysed using generalised 
linear models (GLI^) following the equation STIC = 1 + 
5 genotypes nested within fixed effects + error implemented 
via the S-PLUS software. The same model is used for 
STADH. Fixed effects that were considered were breed, 
finish (Domestic, Korea, Japan), cohort, region (pasture v 
grain, north v south) and the covariate of age. Age was 

10 included since LOX affects cross-linking of collagen and 
cross-linking is an age related process, with cross- 
linking increasing over time. The genotypes were nested 
within region and breed since pure-bred offspring of 
taurine sires were not pastured in the north. The size of 

15 the effects associated with genotype was estimated by the 
comparison of variances (eg, Anders son- Eklund and Rendel, 
1993). To estimate the size of effect associated with 
genotypic substitution, the same model was fitted without 
the LOX genotypes. Residuals were extracted and compared 

2 0 to the LOX genotypes. These were analysed using an 

analysis of variance to obtain adjusted means for each 
genotype . 

Results 

25 There are two alleles (Figure 7) for the LOX 

polymorphism and both these alleles are found in all the 
breeds, although there are clear differences in genotype 
frequency within the breeds. There is no consistent 
difference between zebu and taurine breeds in frequency of 

3 0 the genotypes (Table 12) . The Hereford breed differs 

radically in genotype frequencies to all the other breeds 
in the sample. It has high frequencies of genotype x 22' 
while all other breeds have high frequencies of genotype 

3 5 The STIC and STADH values are correlated with 

R=0.52 (Figures 8 - 10). The plots indicate that while 
the STIC values show two clear extremes the STADH values 
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have only a long tail and do not show two discrete 
distributions. This reflects that the sample was selected 
only on STIC . 

The analyses (Tables 13 and 14) of LOX against 
5 STIC and STADH show consistent statistically significant 
associations. The first and the second samples as well as 
the combined samples of both STIC and STADH show 
associations to LOX genotypes at P<0.05. The association 
of LOX appears stronger to STADH than to STIC. The 

10 association in the second sample of STADH phenotypes has 
extremely high statistical significance but this may be 
due to sampling in small populations and the congi»uence of 
extreme phenotypes with particular genotypes. The 
combined analysis of STADH is less extreme than the second 

15 sample but shows confirmatory linkage to LOX (P<0.01). 

Nevertheless, it is clear that these are not 
large genotypic substitution effects and some analyses do 
not show statistical significance. When the LOX genotypes 
were compared to residual STIC and STADH, none of these 

2 0 associations was statistically significant, whether by 

sample or data combined (Table 15) . Interestingly, some 
of the comparisons show that the heterozygotes are of 
extreme phenotype, opening the possibility of 
overdominance at this locus. This will need to be 

2 5 confirmed using other polymorphisms at the LOX gene that 

show larger genotypic substitution effects. 

Discussion 

Consistent with those earlier analyses, the STADH 

3 0 values show greater association to the DNA marker than the 

STIC values, even though the samples are extreme for STIC, 
with STADH values only more dispersed than normal due to 
the correlation between traits (Figure 9). STIC was 
chosen upon which to select extremes rather than STADH. 
35 However, both of these measurements evaluate aspects of 
connective tissue strength, the adhesion measures the 
force required, in crude terms, to pull a muscle apart, 
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the force applied perpendicular to the fibre bundles, 
while the instron compression measures how much the muscle 
can be flattened without being torn or cut. These are not 
perfectly correlated as can be seen by inspection of the 
5 distribution of STIC and STADH values (Figure 10) . 

Example 5 

The association between the marker and STIC was 
examined in Example 4 using two batches of extreme 
10 animals. The results show that there are significant 

associations between the genotypes of the marker and STIC 
(instrom compression, P < 0.05), and STADH (adhesion, P 
<0.01). The results suggest that the gene LOX either 
contains or is located near the genetic factor associated 
15 with connective tissue strength. 

Because this study was carried out on a relative 
small population (253) with extreme animals only, the same 
marker was tested on different populations to see if the 
association is still valid. 

20 

Materials and Methods 

In addition to the population, there are two 
other groups containing animals chosen from the two tails 
of instron compression (LDIC, 13 6) and peak force (LDPF, 
25 131) for the LOX gene study. These three extreme groups 
together with 559 non-extreme individuals form the base 
for these analyses on the LOX marker. A total of 917 
individuals were used for the study (Table 16) . 

Due to the nature of the populations, the 
3 0 analyses were carried out to the three datasets. 

Extreme animals only (389). The extreme animals 

from LDIC, LDPF and STPF were pooled together. 

Non-extreme animals (559). 

Combined data (917). The combination of 1 and 2. 
3 5 Beside the traits STIC, LDPF and LDIC, a range of 

other traits was also evaluated to see if there is any 
effect of LOX gene on other meat quality traits (Table 
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17) . The intramuscular fat measurements from LD_FAT% and 
NIR_FAT% were combined to make a single trait. 

The mixed model procedure (MLX) in SAS (version 
8.0) was used to .run the statistical analyses. The fixed 
5 effects in the model include finish group and LOX marker. 
Sire and contemporary groups are treated as random 
effects. All these effects were nested within individual 
breeds. The statistical model used is as follow: 

Trait = mean + sire within breed + contemporary 
10 group within breed + finish within breed +LOX within breed 
+Carcass weight. 

Contemporary group was defined as the 
combination of herd of origin, cohort and kill code. The 
individuals without electrical stimulation were removed 
15 from the analysis data. Carcass weight is being used as a 
covariate to adjust for the age difference. 

A full contrast model would be performed if a 
significant marker-trait association was identified from a 
mixed model (or GLM) analysis. The purpose of conducting 
2 0 such the test is to further examine the possibility of 
additive or dominance or overdominance effect among the 
genotypes of the LOX marker. The full contrast of 3 
genotypes (11, 12 and 22) is set up in SAS as follow: 
contrast 'Additive Test' lox(bcode) 1 0 -1; 

2 5 Contrast 'Homozygotel vs Heterozygote ' 1-10 

Contrast 'Heterozygote vs Homozygote2' 0 1 -1 
contrast 'Dominance Test' lox(bcode) -12-1; 
contrast 'Recessive Test' lox(bcode) -1 -1 2; 
contrast ' OverDominance Test' lox(bcode) 2 -1 -1; 

3 0 For the extreme population in which the animals 

with extreme phenotypes were genotyped, multi-trait 
logistic regression method (Henshall and Goddard, 1999) 
was applied to take the potential correlation of traits 
into account. The method is regression based, but instead 
3 5 of regressing phenotype on genotype , the regression is 
genotype on phenotype. This replaces the assumption that 
phenotypes are unselected with the assumption that there 
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was no selection based on genotypes (Henshall and Goddard, 
1998) . Prior to using logistic regression method, MLX 
model was used to all data (917 animals) to derive 
predicted values , of individual animals. The predicted 
5 phenotype values for the extreme animals after adjusting 
for significant fixed effects were then used for Logistic 
regression analyses. The analyses started with single 
trait logistic regression method and then proceeded to 
multi- trait logistic regression method. 

10 The genotype frequency distribution of the 

marker in different populations is shown in Table 18. From 
the table, it can be seen that the Hereford breed differs 
remarkably in genotype frequencies to all the other breeds 
in the populations. In order to investigate the potential 

15 effect of skewed genotypes of Hereford breed on the 
overall results, a set of additional analyses were also 
pursued to the populations by removing the Hereford 
individuals from the data sets. 

Results and Discussions 
20 Part I . Extreme Animals (Table 19) 

Extreme animals for LDIC 

The first test was conducted to the sample 

containing the selected animals for LDIC (136). The 

25 results from the analysis of variance reveal that there 

was no close association between any genotypes of LOX 
marker and LDIC. The same conclusion was held to other 
meat quality traits. 

3 0 Extreme animals for LDPF 

Like LDIC sample, there was no significant 
variation detected between the LOX marker and any meat 
quality trait in the batch animals selected for LDPF. The 
results are not surprising as the initial QTL for 

35 tenderness in CBX experiment was identified in instron 
compression measurement of Semitendinosus muscles. 
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Combined extreme animal data 

Analysis of Variance. As sire effect was 

confounded with other effects, it had to be removed from 
the model and v GLM (generalised linear model) was 
5 performed. In this case, contemporary group was treated as 
a fixed effect rather than a random effect. The analysis 
of variance has shown that out of 21 meat quality traits 
tested, STIC and LDL had significant results (P < 0.05). 

Full Contrast Model. The results from full 

10 contrast model are given below. For STIC, it can be seen 
that there was no additive effect between the two 
homozygous genotypes (11 and 22) . However, the highly 
significant difference between the phenotypes of 11 and 12 
obviously contributed to the detection of dominant and 

15 overdominant effects. Nothing was remarkable for LDL. 

Logistic Regression. After adjusting for the 
significant fixed effects on all data, logistic regression 
was applied to STIC. Multi-trait logistic regression model 
was also fitted to take the potential correlations between 

20 ST measurements into account (STIC, STPF and STADH). The 
results confirm the findings from the other methods. That 
is, LOX genotypes did have a correlation with STIC. The 
regression co-efficiency between lox marker and STIC is 
shown in the output of logistic procedure (below) . The 
' 25 allele substitution effect of the lox marker could be 
derived from the formulae suggested by Henshall and 
Goddard (1999) based on the total variance of whole data. 
The multi-trait logistic regression test on STIC, STPF and 
STADH has shown that both STIC and STPF had significant 

3 0 effects on LOX gene marker. STPF was marginally non- 
significant in GLM analysis. (Table 20) 

Part II. CRC Non-Extreme Animals 

The non-extreme animals (559) were genotyped 
35 against LOX marker in CRC I and but were not part of the 
animals involved in marker evaluation Phases III. The 
mixed model analyses of variance show that beside STIC, 
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the significant marker-trait association was also detected 
to the intramuscular fat (FAT) and LDPH. However, full 
contrast test to STIC and FAT has failed to pinpoint the 
genotype causing the significant results. In the case of 
5 LDPH, it seems that 22 genotype had an important role in 
determining the outcomes. (Table 21) 

Part III. Combined Data 

When extreme and non-extreme animals were pooled 

10 together, the results from mixed model analysis of 

variance show that again the lox marker was associated 
with STIC (P < 0.05). The significant results were also 
found in STL, which is the measurement of darkness of 
carcass meat colour. However in both cases, full contrast 

15 model had not be able to identify the significant genotype 
sources. (Table 22) 

Part IV. Removing Hereford individuals from the combined 
population 

2 0 In order to test the possible effect of skewed 

distribution of lox genotypes of Hereford breed, the 
additional analyses were also performed to the combined 
data with the removal of Hereford breed. The results 
indicate that the removal of Hereford animals has changed 
25 little to the overall significant results of STIC in the 
combined population. From the genotype frequency 
distribution table, it can be seen that the majority of 
Hereford individuals were from the three extreme 
populations except one animal from non-extreme CRC 

3 0 population. (Table 23) 

The overall results from the investigation of 
LOX gene effect on meat quality traits have been 
consistent across three populations (extreme, non-extreme 
and combined) . That is, there is a strong association of 
3 5 LOX gene marker with the instron compression measurement 
of Semitendinosus muscles (P < 0.05). The significant 



WO 02/064820 



- 28 - 



PCT/AU02/00122 



results from other meat quality traits vary from one 
population to another. 

Industrial Applicability 
5 The invention is useful in allowing selection 

and breeding of animals which yield more tender meat. 
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Table 1 



Characteristics of the first Cattle Sample 



10 



15 



20 



25 



Total: 169 



Breeds : 



Regions : 



Markets : 



Cohorts : 



Sires : 



83 high peak force 
86 low peak force 
29 Santa Gertrudis 

25 Hereford 

26 Angus 

27 Belmont Red 
31 Brahman 

31 Shorthorn 

3 8 Pasture South 

28 Pasture North 
57 Grain South 
41 Grain North 
72 Korean 

67 Domestic 
2 5 Japanese 
27 Cohorts 

Median: 5 steers per cohort 
bottom quartile: 2 steers per cohort 
top quartile: 9 steers per cohort 
112 sires 

Median: 1 steer per sire 

bottom quartile: 1 steer per sire 

top quartile: 2 steers per sire 
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Table 2 

Distribution of CAST genotypes in the breeds in the first 
sample • v 

5 

Breed Genotype 



11 12 22 23 

10 Angus 0 7 19 0 

Belmont Red 0 8 19 0 

Brahman 6 13 10 2 

Hereford 0 5 17 0 

Santa Gertrudis 3 5 19 0 

15 Shorthorn 0 4 23 0 
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Table 3 

Associations between calpastatin genotypes (castg) and 
tenderness . 

5 

A. Calpastatin by itself 



Analysis of Deviance Table 

Gaussian model 

Response: peakforce 

15 Terms added sequentially (first to last) 

Df Deviance Resid. Df Resid. Dev F Value Pr(F) 

NULL 155 864.6307 

20 castg 3 70.89899 152 793.7317 4.52573 0.004536025 



25 

B. Breed x Calpastatin Interactions 



Analysis of Deviance Table 

30 

Gaussian model 
Response: peakforce 
35 Terms added sequentially (first to last) 



40 





Df 


Deviance 


Resid. Df 


Resid . 


. Dev 


F 


Value 




Pr ( F) 


NULL 






155 


864 . 


.6307 










finish 


2 


101.8455 


153 


762 . 


.7852 


15. 


.97637 


0. 


.0000008 


cohort 


25 


273 . 5882 


128 


489. 


.1971 


3 . 


.43339 


0, 


.0000041 


region 


3 


59.6620 


125 


429. 


.5350 


6. 


.23940 


0. 


.0005889 


breed 


4 


10.8711 


121 


418 


.6639 


0. 


.85267 


0 


.4948672 


castg 


3 


28.2709 


118 


390, 


.3930 


2. 


.95654 


0, 


.0355308 


breed : castg 


6 


33 . 4066 


112 


356, 


.9864 


1 


.74682 


0 


.1166518 



45 
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Table 3 (continued) 

C. Calpastatin genotypes nested with breed and region 

Analysis of Deviance Table 
Gaussian model 
Response: peakforce 

Terms added sequentially (first to last) 

15 Df Deviance Resid. Df Resid. Dev F Value 

Pr { F) 

NULL 155 864.6307 

finish 2 101.8455 153 762.7852 18.28842 

0.0000002 

20 cohort 25 273.5882 128 489.1971 3.93027 

0.0000005 

region 3 59.6620 125 429.5350 7.14235 

0.0002128 

breed in region 7 26.8943 118 402.6408 1.37983 

25 0.2219802 

castg in (region/breed) 17 121.4139 101' 281.2269 2.56498 

0.0018938 



30 
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Table 4 



Analysis of CAST against residual peakforce measurements (XI) . 

5 Analysis of Deviance Table 
Gaussian model 



Response: XI 

10 

Terms added sequentially (first to last) 

Df Deviance Resid. Df Resid. Dev F Value Pr(F) 

15 NULL 154 402.3160 

castg 3 21.97947 151 380.3365 2.90874 0.03654659 



Call: 

20 

glm(formula = XI castg, data = calppftest, na. action = na.omit) 
Coefficients : 



25 (Intercept) castgl castg2 castg3 

0.07693593 -0.4095948 -0.3102917 -0.3504961 



30 Degrees of Freedom: 155 Total; 151 Residual 



Residual Deviance: 380.3365 



35 

Model from which residuals were calculated 



Analysis of Deviance Table 
40 Gaussian model 



Response: peakforce 

Terms added sequentially (first to last) 

45 





Df 


Deviance Resid, 


. Df 


Resid. 


. Dev 


F 


Value 




Pr (F) 


NULL 








160 


906. 


.5175 










finish 


2 


106 


.0217 


158 


800. 


.4958 


15. 


.52293 


0 . 


.0000010 


cohort 


26 


290 


.8387 


132 


509. 


.6571 


3 . 


.27558 


0 . 


.0000057 


f inlwt 


1 


8 


.0662 


131 


501 . 


.5908 


2 . 


.36200 


0. 


.1269337 


region 


3 


63 


.9431 


128 


437. 


.6478 


6 


.24139 


0 


.0005629 


region 


7 


24 


.4323 


121 


413. 


.2154 


1. 


.02206 


0 


.4192678 



55 

glm( formula: peakforce = finish + cohort + f inlwt + region/breed, 



data = calppf, na . action = na.omit) 
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Table 5 

Analysis of Variance tables between CAST genotypes and residual 
peak force measures (XI) along with the table of adjusted means 
5 associated with each genotype. 

Analysis of Variance Table 

Response: XI 

10 

Terms added sequentially (first to last) 

Df Sum of Sq Mean Sq F Value Pr(F) 
15 castg 3 21.9795 7.326490 2.90874 0.03654659 

Residuals 151 380.3365 2.518785 



20 



30 



Tables of adjusted means 



25 Grand mean 

0.076936 
se 0.318703 



castg 



ell cl2 c22 c23 

1.1473 0.3281 -0.1932 -0.9746 
35 se 0.5290 0.2479 0.1564 1.1222 
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Table 6 

Characteristics of the second sample of 77 animals. 





r Pr* t"al • 
1 Otai . 


in 
/ / 


V 


-J 






mgn pea.K lorce 






3 o 


low peak force 




Breeds : 


11 


Belmont Reel 






11 


Hereford 






13 


Brahman 


1U 




13 


Shorthorn 






14 


Santa Gertrudis 






15 


Angus 




Regions : 


24 


Pasture South 






12 


Pasture North 


15 




21 


Grain South 






20 


Grain North 




Markets : 


35 


Korean 






25 


Domestic 






17 


Japanese 


20 


Cohorts : 


22 


Cohorts 



Median: 3 steers per cohort 
bottom quartile: 2 steers per cohort 
top quartile: 5 steers per cohort 
Sires: 64 sires 

25 Median: 1 animal per sire 

bottom quartile: 1 animal per sire 
top quartile: 1 animal per sire 
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Table 7 

Distribution of CAST genotypes in the second sample 



5 



Breed Genotype 



11 12 22 

Angus 0 3 12 

Belmont Red 0 3 8 

10 Brahman 3 7 3 

Hereford 0 3 8 

Santa Gertrudis 1 4 9 

Shorthorn 0 0 13 
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Table 8 



Associations between calpastatin genotypes and LD peak force in 

i 

the second sample. 

5 Analysis of Deviance Table 
Gaussian model 
Response: ldpeakforce 

10 

Terms added sequentially (first to last) 

Df Deviance Resid. Df Resid. Dev F Value Pr(F) 

15 NULL 76 205.9332 

castg 2 17.90313 74 188.0300 3.522925 0.03455689 

' Coefficients : 

(Intercept) castgl castg2 
20 5.227591 -0.1205 -0.3719088 



Degrees of Freedom: 77 Total; 74 Residual 
2 5 Residual Deviance: 188.03 



3 0 Single term deletions 
Model : 

ldpeakforce = Islortwait + buttemp + finish + cohort + region + 
3 5 breed + castg 



Final Call: 

40 glm( formula = ldpeakforce castg, data = calppfr, na. action = 
na . omit ) 



Coefficients : 
45 (Intercept) castgl castg2 

5.195064 -0.07055556 -0.37438 



Degrees of Freedom: 67 Total; 64 Residual 
50 Residual Deviance: 161.5878 
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Table 9 

The amount of each haplotype found between the alleles of the 
casts microsatellite and the cast3 D/E SNP on both cattle samples. 

Twenty- six haplptypes were found in animals that are homozygous 
5 for one or the other locus . 

allele 

haplotype cast3 cast5 amount 



10 



15 



20 



25 



30 



35 



1 


al 


ml 


3 


2 


al 


zn2 


0 


3 


al 


m3 


14 


4 


al 


m4 


0 


5 


al 


m5 


0 


6 


al 


m6 


3 


7 


al 


m7 


6 


8 


al 


m8 


2 


9 


al 


m9 


6 


10 


al 


mlO 


5 


11 


al 


mil 


0 


12 


al 


ml2 


1 


13 


al 


ml 3 


0 


14 


a2 


ml 


0 


15 


a2 


m2 


1 


16 


a2 


m3 


108 


17 


a2 


m4 


1 


18 


a2 


m5 


4 


19 


a2 


m6 


3 


20 


a2 


m7 


42 


21 


a2 


m8 


2 


22 


a2 


m9 


110 


23 


a2 


mlO 


17 


24 


a2 


mil 


6 


25 


a2 


ml 2 


1 


26 


a2 


ml 3 


1 
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Table 10 

A heterogeneity test for associations between alleles at CAST5 
with alleles at CAST3 D/E. 

5 Analysis of Deviance Table 



Poisson model 



Response: score 



10 

Terms added sequentially (first to last) 

Df Deviance Resid. Df Resid. Dev Pr(Chi) 

NULL 25 926.9836 

cast3 1 220.4994 24 706.4842 0.0000000000 

15 casts .12 671.7497 12 34.7345 0.0000000000 

cast3:cast5 12 34.7344 0 0.0001 0.0005161421 



Table 11 

20 Tests for association between CAST 5 microsatellite and LD peak 

force measurements in both cattle samples. 

Part A. Calpastatin by itself 

25 Analysis of Deviance Table 

Gaussian model 

Response: ldpf 

30 

Terms added sequentially (first to last) 

Df Deviance Resid. Df Resid. Dev F Value Pr(F) 
NULL 491 2266.640 

malll 12 136.0104 479 2130.629 2.548113 0.002843108 

35 

Part B. Breed by calpastatin interactions 



Analysis of Deviance Table 



WO 02/064820 



PCT/AU02/00122 



- 40 - 

Table 10 (Continued) 
Gaussian model 
5 Response: Id 



Terms added sequentially (first to last) 



10 



15 





Df 


Deviance 


Resid. Df 


Resid . 


Dev 


F 


Value 




Pr(F) 


NULL 






491 


2266 


.640 










market 


2 


218.4756 


489 


2048 


.164 


42 


.35776 


0 


.0000000 


cohort 


26 


706.3980 


463 


1341 


.766 


10 


.53504 


0 


.0000000 


finish 


3 


97 .7159 


460 


1244 


.050 


12 


. 63003 


0 


.0000001 


breed 


4 


23.5779 


456 


1220 


.472 


2 


.28562 


0 


.0594827 


cast5 


12 


65.7617 


444 


1154 


.711 


2 


.12497 


0 


.0146070 


breed : cast 5 


26 


76.7170 


418 


1077 


.994 


1 


. 14414 


0 


.2865614 



Part C. Calpastatin genotypes nested with breed and finish (region) 
Analysis of Deviance Table v 

20 ■ 

Gaussian model 



Response: Id 

2 5 Terms added sequentially (first to la 

Df Deviance Re 

NULL 

market 2 218.4756 
cohort 26 706.3980 

3 0 breed 4 30.7280 

finish %±n% reed 6 111.2067 

casts %±n% (breed/ finish) 65 211.8811 



id. Df Resid. Dev F Value Pr(F) 
491 2266.640 

489 2048.164 42.90121 0.00000000 
463 1341.766 10.67020 0.00000000 
459 1311.038 3.01697 0.01799825 
453 1199.832 7.27908 0.00000022 
388 987.950 1.28019 0.08307834 
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Table 11 

Characteristics of the Cattle Sample 



Total: 



high instron compression 
low instron compression 



166 

87 

89 



87 
39 
38 



Breeds : 

Angus 

10 Belmont Red 

Brahman 
Hereford 
Santa Gertrudis 
Shorthorn 

15 

Regions : 

Pasture South 
Pasture North 
Grain South 
20 Grain North 



25 
25 
33 
32 
26 
25 



47 

39 
43 
37 



12 
12 
18 
10 
15 
17 



20 
21 
20 
26 



25 



Markets : 

Korean 

Domestic 

Japanese 



81 
47 
38 



22 
45 
22 



30 



Cohorts: 2 5 14 

Median: steers per cohort 6 3 

bottom quartile: 3 1 

top quartile: 10 11 



Sires : 



35 



Median: steers per sire 
bottom quartile: 
top quartile: * 



113 

1 

1 

2 



62 
1 
1 
2 
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Table 12 

Distribution of LOX genotypes in the breeds in the combined 
sample . 

5 Breed " Genotype 

11 12 22 

Angus 12 16 5 

Belmont Red 19 14 2 

Brahman 20 21 4 

10 Hereford 1 7 27 

Santa Gertrudis 18 5 1 

Shorthorn 23 11 3 
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Table 13 

Associations between LOX genotypes (loxg) and STIC . 
A. First Sample 

5 

Analysis of Deviance Table 
Gaussian model 
10 Response: stic 

Terms added sequentially (first to last) 



Df Deviance Resid.Df Resid. Dev F Value Pr(F) 



NULL 








144 


66. 


41782 










market 




2 


4.68851 


142 


61. 


72932 


10 


.11346 


0 


.0001124 


age 




1 


0.33061 


141 


61. 


39871 


1 


.42631 


0 


.2356143 


cohort 




23 


25.50396 


118 


35. 


89474 


4 


.78382 


0 


.0000000 


region 




3 


2,05111 


115 


33. 


84363 


2 


.94960 


0 


.0371506 


breed %in% 


region 


8 


3.73686 . 


107 


30. 


10677 


2 


.01517 


0 


.0537829 


loxg %in% 


(region/breed) 20 


9.94057 


87 


20. 


16620 


2 


.14425 


0 


.0081786 



B. Second Sample 



25 Analysis of Deviance Table 



Gaussian model 
Response: stic 

Terms added sequentially (first to last) 
3 0 Df Deviance Resid.Df Resid. Dev F Value Pr(F) 



NULL 










77 


36 


.73118 










market 




2 


18. 


56015 


75 


18 


.17103 


50. 


57448 


0 


.0000000 


age 




1 


1. 


22303 


74 


16 


.94800 


6. 


66525 


0 


.0131586 


region 




3 


0. 


49228 


71 


16 


.45571 


0. 


89428 


0 


.4514796 


breed %in% 


region 


9 


2. 


26136 


62 


14 


.19435 


1. 


36933 


0 


.2303598 


loxg %in% 


(region/breed) 17 


5. 


93715 


45 


8 


.25720 


1. 


90331 


0 


.0433708 
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NOTE: cohort could not be fitted as it required a model with more 
terms than degrees of freedom. cohort was dropped since that allowed 
a maximum of other terms to be fitted. 

•» 

5 C. Combined Sample 

Analysis of Deviance Table 



Gaussian model 
10 Response:, stic 

Terms added sequentially (first to last) 

Df Deviance Resid.Df Resid. Dev F Value Pr(P) 

NULL 222 105.0131 

market 2 9.28932 220 95.7238 14.63402 0.0000015 

15 age 1 0.66404 219 95.0598 2 . 0922 1 . 0 . 150007 9 

cohort 23 28.46167 196 66.5981 3.89890 0.0000002 

region 3 0.97412 193 65.6240 1.02306 0.3840448 

breed %in% region 9 2.08684 184 63.5372 0.73056 0.6804229 

loxg %±tx% (region/breed) 24 12.75507 160 50.7821 1.67448 0.0327546 

20 
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Table 14 



Associations between LOX genotypes (loxg) and STADH. 
5 A. First Sample 



Analysis of Deviance Table 
Gaussian model 
Response: stadh 
10 Terms added sequentially (first to last) 

Df Deviance Resid.Df Resid. Dev F Value Pr(F) 



NULL 










139 


4. 


118160 










market 




2 


0 


.045611 


137 


4. 


072549 


1 


.63665 


0. 


.2008490 


age 




1 


1 


.030876 


136 


3 . 


041674 


73 


.98193 


0. 


.0000000 


cohort 




23 


1 


.131428 


113 


1. 


910246 


3 


.53035 


0. 


.0000129 


region 




3 


0 


.053546 


110 


1 . 


856699 


1 


.28094 


0, 


,2863452 


breed %±n% 


region 


8 


0 


.192061 


102 


1. 


664639 


1 


.72293 


0. 


.1050769 


loxg %in% 


(region/breed) 19 


0 


.508104 


83 


1. 


156535 


1 


.91919 


0. 


.0229812 



20 

B . Second Sample 



Analysis of Deviance Table 



2 5 Gaussian model 
Response: stadh 

Terms added sequentially (first to last) 



Df Deviance Resid.Df Resid. Dev F Value Pr(F) 



NULL 








76 


3.362345 










market 




2 0 


.560443 


74 


2.801903 


25 


.18352 


0 


.0000000513 


age 




1 0 


.486387 


73 


2.315516 


43 


.71164 


0 


.0000000425 


region 




3 0 


.390962 


70 


1.924554 


11 


.71193 


0 


.0000090738 


breed %in% 


region 


9 0 


.413632 


61 


1.510922 


4 


.13035 


0 


.0006637134 


loxg %xn% 


(region/breed) 




















17 1 


.021326 


44 


0.489595 


5 


.39922 


0 


.0000032684 



WO 02/064820 



- 46 - 



PCT/AU02/00122 



C. Combined Sample 

Analysis of Deviance Table 

•. 

5 Gaussian model 
Response: stadh 

Terms added sequentially (first 

Df Deviance 

10 NULL 

market 2 0.2 68508 

age 1 1.252352 

cohort 23 2.111912 

region 3 0.065737 

15 breed %±ti% region 9 0.325818 
loxg %±n% (region/breed) 

23 0.784128 



to last) 

Resid.Df Resid. Dev F Value Pr(F) 



216 


7 


.480742 










214 


7 


.212234 


7 


.78709 


0 


.0005990 


213 


5 


.959882 


72 


.63985 


0 


.0000000 


190 


3 


.847971 


5 


.32594 


0 


.0000000 


187 


3 


.782233 


1 


.27098 


0 


.2863505 


17 8 


3 


.456415 


2 


.09982 


0 


.0326124 


155 


2 


.672287 


1 


.97746 


0 


.0079540 
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Table 15 

Estimated sizes of effects of genotype substitutions at LOX on 
instron compression and adhesion of the semi tendinosus muscle. 



10 



15 



Response: stic. resid 

Grand-mean 0.027133 se 

Loxg 111 112 

0.03958 -0.02321 

se 0.07654 0.07897 

Response: sticbox3 . resid 

Grand-mean -0.012888 se 

Loxg 111 112 

0.033317 -0.058769 

se 0.055916 0.058431 



0.053648 
122 
0.06503 
0.11751 



0.039137 
122 

-0.013213 
0.085115 



20 



Response: sticfull . resid 

Grand-mean -0.012888 se 0.039137 

loxg Hi H2 122 

0.033317 -0.058769 -0.013213 

se 0.055916 0.058431 0.085115 



Response: sticadh. resid 

Grand-mean 0.002521 se 0.010321 

25 Loxg 111 112 122 

-0.001426 -0.002304 0.011292 

se 0.014860 0.015390 0.022384 



Response: sticbox3adh . resid 

3 0 Grand-mean 0.004308 se 0.012879 

loxg Hi H2 122 

-0.003154 -0.007446 0.023524 

se 0.018594 0.018891 0.028111 



35 
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Response: st icf ulladh. resid 

Grand-mean 0.0037399 se 0.0092835 

loxg in 112 122 

-0.009303 -0.005244 0.025767 

5 se 0.013379 0.013763 0.020180 
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Table 16 



Information on the data sets 



Effect 


Class 


LDIC 


LDPF 


Combined 
Extreme 


Non- 
extreme 


Combined 


Total 




136 


131 


398 


543 


916 


Sires 




96 


96 


171 


61 


227 


Cohorts 




24 


24 


26 


10 


30 


Breeds 


Angus 


19 


22 


62 


134 


196 




Belmont Red 


23 


25 


66 


140 


200 




Brahman 


24 


27 


76 


73 


142 




Hereford 


27 


14 


67 


1 


68 




Santa Gertrudis 


25 


24 


73 


195 


257 




Shorthorn 


18 


19 


53 


0 


53 


Regions 


Pasture South 


27 


27 


92 


93 


185 




Pasture North 


24 


25 


78 


144 


212 




Grain South 


54 


41 


125 


168 


291 




Grain North 


31 


38 


102 


138 


228 


Markets 


Domestic 


47 


51 


130 


240 


361 




Korean 


61 


58 


189 


201 


376 




Japaness 


28 


22 


78 


102 


179 
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Table 17 



Meat quality traits tested for LOX gene marker 



Code 


Tra.it - 




Intramuscular Fat percentage (Soxhylet Method) 


LD IC 


Longissimus dorsi Instrom compression 




Longissimus dorsi initial yield (Nth kills only) 




Longissimus dorsi cooking losses 




Longissimus dorsi Peak Force -must use "Stim" also 


T OP-TV 


Longissimus dorsi Peak Force - initial yield (Nth) 


T Tl -=» 


Longissimus dorsi a* colour 


T n Vk 

JbJJ D 


Longissimus dorsi b* colour 


T T\ 1 
LiU A 


Longissimus dorsi L* colour 


LD pH 


Longissimus dorsi ultimate pH 


NXK. r at-s 


Intramuscular Fat percentage (NIR method) 


ST AdhRS 


oomi.LciluxiiObUS oxior tnose aQJleSlOIl 


ST_IC 


Semitendinosus Instrom compression 


ST_IY 


Semitendinosus initial yield (Nth kills only) 


ST_LOSS 


Semitendinosus cooking loss% 


ST_PF 


Semitendinosus Peak Force 


ST_PF-IY 


Semitendinosus Peak Force - initial yield (Nth) 


ST_a 


Semitendinosus a* colour 


ST_b 


Semitendinosus b* colour 


ST_1 


Semitendinosus L* colour 


ST_pH 


Semitendinosus ultimate pH 


TenderQ 


Tenderness Quality as measured by PF (x 100) 



5 
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Table 18 

Distribution of LOX genotypes in the breeds in the three datasets 



Breed 


Extreme 


Non- extreme C 


Combined 




11 


12 


22 


11 


12 


22 


11 


12 


22 


Angus 


21 


33 


7 


38 


75 


21 


59 


109 


28 


Brahman 


35 


35 


6 


19 


40 


14 


52 


71 


19 


Belmont Red 


31 


23 


12 


59 


60 


21 


87 


83 


30 


Hereford 


4 


17 


46 


0 


1 


0 


4 


18 


46 


Santa Gertrudis 


42 


29 


2 


120 


62 


13 


156 


86 


15 


Shorthorn 


34 


18 


1 


0 


0 


0 


34 


18 


1 


Total 


167 


156 


74 


236 


238 


69 


392 


385 


139 



5 
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Table 19 



The GUI Procedure 



10 



15 



Dependent Variable: STIC 



<.0O01 



Model 

1 

Error 

Corrected Total 



DF 



297 



Sum of 
Squares 



90.3108958 



97 11.5922526 
394 101.9031484 



Mean Square 



0.3040771 



0.1195078 



P Value 



2.54 



20 



R-Square 
0.886242 



Coeff Var 
16.41081 



Root MSE 
0.345699 



STIC Mean 
2.106532 



25 



Source 



DF 



Type III SS 



Mean Square 



F Value 



Pr 



con temp {Be ode) 



30 



35 



<.0001 



Fingp(Bcode) 



0.2004 



0.1199 



lox ( Bcode ) 



0.0145 



0.0243 



269 
3 
1 
12 
1 



72.05453149 
0.56463473 
0.29426415 
3.23498136 
0.62547588 



0.26786071 
0. 18821158 
0.29426415 
0.26958178 
0.62547588 



2.24 
1.57 
2.46 
2.26 
5.23 



40 



Contrast 



DF 



Contrast SS 



Mean Square 



F Value 



Additive Test 11-22 



45 



50 



55 



60 



65 



70 



0.0636 



11-12 



<.0001 



12-22 



0.2283 



Dominance Test 



0.0052 



Recessive Test 



0.7770 



OverDomi nance Test 



0.0004 



Parameter 

Additive Test 11-22 

11- 12 

12- 22 

Dominance Test 
Recessive Test 
Over Dominance Test 



LDL 



Dependent Variable: LD1 



0.42074192 
2 .17031475 
0.17567132 
0.97536352 
0.00964320 
1.58776355 



Estimate 

-0.54698783 
-0.93226650 
0.38527867 
1.31754517 
0.16170916 
-1.47925433 



0.42074192 
2. 17031475 
0. 17567132 
0.97S36352 
0.00964320 
1.58776355 



Standard 
Error 

0.29151964 
0.21876446 
0.31777709 
0.46119045 
0.56927504 
0.40583359 



3.52 
18.16 
1.47 
8.16 
0.08 
13.29 



t Value 

-1.88 
-4.26 
1.21 
2.86 
0.28 
-3.64 



Pr > |t| 

0.0636 
<.00Ol 
0.2283 
0.0052 
0.7770 
0.0004 



The GLM Procedure 
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Table 19 (Continued) 



Source 



Sum of 
Squares 



Mean Square F Value Pr > 



Model 



295 



4070.356476 



13.797819 



1.96 



10 



Error 

Corrected Total 



96 
391 



675.059334 
4745.415310 



7.031873 



15 



R-Square 
0.857745 



Coeff Var 
6.920726 



Root MSE 
2.651768 



LD1 Mean 
38.31633 



20 



Source 



DF 



Type III SS 



Mean Square 



F Value 



25 



30 



0.0278 



contemp ( Be ode ) 

3 

Fingp(Bcode) 



0.3161 



Stim 



0.8678 



lox(Bcode) 



0.0427 



0.3590 



267 
3 
1 
12 
1 



2628.383772 
25.195055 
0. 195851 • 
160.954945 
5.974613 



9 .844134 
8 .398352 
0.195851 
13.412912 
5.974613 



1.40 
1.19 
0. 03 
1.91 
0.85 



35 



Contrast 



DF 



Contrast SS 



Mean Square 



F Value 



Pr > 



Additive Test 



40 



45 



50 



55 



11 vs 12 



0. 5091 



12 vs 22 



0.7024 



Dominance Test 



0.9600 



Recessive Test 



0.4966 



Over Dominance Test 



0.3129 



Parameter 

Additive Test 11-22 

11- 12 

12- 22 

Dominance Test 
Recessive Test 
Over-Dominance Test 



5.88718974 
3 .08780445 
1.03256209 
0.01778689 
3.27515251 
7 .23682011 



Estimate 

-2.04608559 
-1.11200497 
-0.93408062 
0.17792435 
2.98016622 
-3.15809056 



5. 88718974 
3. 08780445 
1.03256209 
0.01778689 
3.27515251 
7.23682011 



Standard 
Error 

2.23617248 
1.67809812 
2.43759636 
3.53769859 
4.36676923 
3.11305081 



t Value 

-0.91 
-0.66 
-0.38 
0.05 
0.68 
-1.01 



0. 84 
0.44 
0. 15 
0.00 
0.47 
1.03 



Pr > |t| 

0.3625 
0.5091 
0. 7024 
0. 9600 
0.4966 
0.3129 
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Table 20 (continued) 

Single Trait Logistic Regression 



10 



15 



20 



25 



30 



35 



40 



45 



The LOGISTIC Procedure ■ 
Model Information 



Data Sec 

Response Variable 
Number of Response Levels 
Number of Observations 
Link Function 
Optimization Technique 



WORK . EXTRERES 
lox 

3 

389 
Log it 

Fisher's scoring 



Ordered 
Value 

1 
2 
3 



Response Profile 



lox 



11 
12 
22 



Total 
Frequency 

165 
153 
71 



NOTE: 8 observations were deleted due to missing values for the response or explanatory 
variables. 

The LOGISTIC Procedure 
Testing Global Null Hypothesis: BETA=0 



Test 


Chi-Square 


DF Pr 


> ChiSq 


Likelihood Ratio 5.0720 


1 


0.024 3 


Score 


5.0502 


1 


0.0246 


Wald 


5.0083 


1 


0.0252 




Analysis of Maximum Likelihood Estimates 






Standard 






Parameter 


DF Estimate Error 


Chi-Square 


Pr > ChiSq 


Intercept 


1 0.9245 0.5608 


2 .7174 


0.0993 


Intercept2 


1 2.7482 0.5777 


22 .6278 


<.0001 


sticpred 


1 -0.5827 0.2604 


5.0083 


0.0252 



50 



Odds Ratio Estimates 



Effect 
sticpred 



Point 
Estimate 

0.558 



95% Wald 
Confidence Limits 



0.335 



0.930 
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Table 20 (continued) 



Multi-trait Logistic Regression 



10 



15 



20 



25 



The LOGISTIC Procedure 
Testing Global Null Hypothesis: BETA=0 
Test Chi-Square DF Pr > ChiSq 



Likelihood Ratio 

Score 

Wald 



Parameter 

Intercept 11 
Intercept 12 
sticpred 
stadhpred 
stpfpred 



14.7234 
14.5359 
13.9424 



0.0021 
0.0022 
0.0030 



Analysis of Maximum Likelihood Estimates 



DF 



Estimate 

-0.9210 
0.9467 
-0.8681 
-0.4025 
0.5689 



Standard 
Error 

0.8702 
0.8711 
0.3191 
0.7724 
0.1936 



Wald 
Chi-Square 

1.1202 
1.1813 
7.4014 
0.2715 
8.6313 



Pr > ChiSq 



2899 
2771 
0065 
6023 
0033 



Odds Ratio Estimates 

3 0 Point 95% Wald 

Effect Estimate Confidence Limits 

sticpred 0.420 0.225 0;785 

stadhpred 0.669 0.147 3.039 

35 stpfpred 1.766 1:208 2.581 
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Table 21 



10 



15 



20 



25 



30 



35 



40 



45 



50 



- 56 - 



The Mixed Procedure 



Model Information 



Data Set 

Dependent Variable 
Covariance Structure 
Estimation Method 
Residual Variance Method 
Fixed Effects SE Method 



WORK.CRC1 
STIC 

Variance Components 

REML 

Profile 

Model-Based 



Degrees of Freedom Method Containment 



Covariance Parameter Estimates 



Standard 



Cov Parm 


Estimate 


Error 


Value 




Pr 2 






SirelD(Bcode) 


0.002544 


0.002147 


1 


18 


0 


1181 






contemp( Bcode) 


0.008813 


0.003896 


2 


26 


0 


0118 






Residual 


0.07009 


0.005174 


13 


55 


< 


0001 








Type 3 Tests 


of Fixed Effects 














Num 


Den 














Effect 


DF 


DF F 


Value 


Pr > F 








Fingp(Bcode) 9 


317 


17.50 


< 


0001 








Stim 


1 


317 


1.75 


0 


1870 








lox(Bcode) 


8 


317 


2.36 


0 


0176 








wt 


1 


317 


6.45 


0 


0115 










Estimates 


















Standard 














Label 


Estimate 


Error 


DF 


t 


Value 


Pr 


> |t| 


Additive Test 11-22 


0.05547 


0.07813 


317 




0.71 




0 


.4782 


11-12 


0. 1044 


0.05660 


317 




1.84 




0 


.0661 


12-22 


-0.04890 


0.07013 


317 




-0.70 




0 


.4861 


Dominance Test 


-0.1533 


0.1007 


317 




-1. 52 




0 


. 1289 


Recessive Test 


-0.00657 


0.1373 


317 




-0.05 




0 


.9619 


OverDominance Test 


0.1598 


0.1170 


317 




1.37 




0 


.1730 



Contrasts 



55 



60 





Num 


Den 






Label 


DF 


DF 


F Value 


Pr > F 


Additive Test 11-22 




317 


0.50 


0.4782 


11-12 




317 


3.40 


0.0661 


12-22 




317 


0.49 


0.4861 


Dominance Test 




317 


2.32 


0.1289 


Recessive Test 




317 


0.00 


0.9619 


OverDominance Test 




317 


1.87 


0.1730 



Instramuscuar Fat 



65 



70 



The Mixed Procedure 
Model Information 



Data Set 

Dependent Variable 
Covariance Structure 



WORK.CRC1 
Fat 

Variance Components 
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Table 21 (continued) 



10 



15 



20 



25 



30 



35 



40 



45 



50 



55 



60 



65 



LDPH 



Estimation 


Method REML 










Residual Variance Method Profile 








Fixed Effects SE Method Model-Based 








Degrees of 


Freedom Method Containment 










Covariance Parameter Estimates 










Standard 




Z 






Cov Parm 


Estimate Error 


Value 




Pr Z 


SirelD(Bcode) 


0.04848 0.04510 


1 


.07 


0. 


1412 


contemp (Bcode) 


0.3111 0.1060 


2 


.93 


0. 


0017 


Residual 


1.2479 0.09791 


12 


.75 


< . 


0001 




Type 3 Tests of Fixed Effects 










Num Den 










Effect 


DF DF F 


Value 


Pr > F 




Fingp (Bcode) 9 306 


7.08 


< 


0001 




Stint 


1 306 


0.31 


0 


5790 




■ lox (Bcode) 8 306 


2.00 


0 


0461 




wt 


1 306 


74 .54 


< 


0001 






Estimates 












Standard 










Label 


Estimate Error 


DF 


t 


Value 


Pr 


Additive Test 11-22 


0.3765 0.3427 


306 




1.10 




11-12 


-0.1207 0.2505 


306 




-0.48 




12-22 


0.4972 0.3029 


306 




1.64 




Dominance Test 


0.6179 0.4377 


306 




1.41 




Recessive Test 


-0.8737 0.5964 


306 




-1.47 




Over Dominance Test 


0.2558 0.5184 


306 




0.49 






Contrasts 












Num Den 










Label 


DF DF 


F Value 


Pr > 


F 


Additive Test 11-22 1 306 


1. 


21 


0.2728 


11-12 


1 306 


0. 


23 


0. 6304 


12-22 


1 306 


2. 


69 


0. 1018 


Dominance Test 1 3 06 


1 . 


99 


0. 1591 


Recessive Test 1 306 


2 . 


15 


0.1439 


Over Dominance Test 1 306 


0. 


24 


0.6220 




The Mixed Procedure 












Model Information 











> Itl 

0.2728 
0.6304 
0.1018 
0.1591 
0.1439 
0.6220 



Data Set 

Dependent Variable 
Covariance Structure 
Estimation Method 
Residual Variance Method 
Fixed Effects SE Method 
Degrees of Freedom Method 



WORK.CRC1 
LDpH 

Variance Components 

REML 

Profile 

Model-Based 

Containment 



70 



Covariance Parameter Estimates 

Standard Z 
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Table 21 (continued) 

Cbv Parm 

con temp (Be ode) 
Residual 



Estimate 

0.002771 
0.007878 



Error 

0.000702 
0.000572 



Value 

3.95 
13.77 



Pr 2 

<.0001 
<.0001 



10 



15 



Type 


3 Tests 


of Fixed 


Effects 






Num 


Den 






Effect 


DF 


DF 


F Value 


' Pr > F 


Fingp (Bcode) 


9 


363 


3 .44 


0.0004 


Stim 


1 


363 


1.14 


0.2854 


1 ox (Bcode) 


8 


363 


3.01 " 


0.0027 


wt 


1 


363 


21.61 


<.0001 



Estimates 







Standard 








Label 


Estimate 


Error 


DF 


t Value- 


Pr > |t| 


Additive Test 11-22 


-0.1089 


0.02648 


363 


-4.11 


<.0001 


11-12 


-0.01502 


0.01922 


363 


-0.78 


0.4351 


12-22 


-0.09393 


0.02376 


363 


-3.95 


<.000l 


Dominance Test 


-0.07891 


0.03416 


363 


-2.31 


0.0215 


Recessive Test 


0.2029 


0.04650 


363 


4.36 


<.0001 


Over Dominance Test 


-0.1240 


0,03971 


363 


-3.12 


0.0019 



The Mixed Procedure 
Contrasts 





Num 


Den 






Label 


DF 


DF 


F Value 


Pr > F 


Additive Test 11-22 




363 


16.93 


<.0001 


11-12 




363 


0.61 


0.4351 


12-22 




363 


• 15.63 


<.0001 


Dominance Test 




363 


5.34 


0.0215 


Recessive Test 




363 


19.04 


<.0001 


OverDominance Test 




363 


9.75 


0.0019 
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Table 22 

STIC 



The Mixed Procedure 
Covariance Parameter Estimates 



10 



15 



20 



25 







Standard 


Z 






Cov Parm 


Estimate 


Error Value 




Pr Z 


SirelD(Bcode) 


0.006141 


0.004815 1 


28 


0 


1011 


contempt Bcode) 


0.05433 


0.01003 5 


42 


< 


0001 


Residual 


0.09190 


0.006489 14 


16 


< 


0001 


Type 3 Tests 


of Fixed Effects 










Num 


Den 








Effect 


DF 


DF F Value 


Pr > F 




Fingp (Bcode) 


11 


369 9.24 




0001 




Stim 


2 


369 4.06 


0. 


0181 




lox (Bcode) 


12 


369 1.90 


0. 


0336 




wt 


1 


369 9.93 


0. 


0018 





Estimates 



30 



35 



40 



Label 

Additive Test 11-22 

11- 12 

12- 22 

Dominance Test 
Recessive Test 
Over Dominance Test 





Standard 






Estimate 


Error 


DF 


t Value 


0.06327 


0.08142 


369 


0.78 


0.07985 


0.05752 


369 


1.39 


-0.01658 


0.07440 


369 


-0.22 


-0.09643 


0. 1052 


369 


-0.92 


-0.04669 


0.1450 


369 


-0.32 


0.1431 


0.1198 


369 


1.20 




Contrasts 







Pr > |t| 

0.4376 
0.1659 
0.8238 
0.3598 
0.7476 
0.2328 



Num Den 

Label DF DF F Value Pr > F 

45 



Additive Test 11-22 


1 369' 


0 


60 


0 


4376 


11-12 


1 369 


1 


93 


0 


1659 


12*22 


1 369 


0 


05 


0 


8238 


Dominance Test 


1 369 


0 


84 


0 


3598 


Recessive Test 


1 369 


0 


10 


0 


7476 


OverDominance Test 


1 369 


1 


43 


0 


2328 



STL 

55 Covariance Parameter Estimates 

Standard Z 





Cov Parm Estimate 


Error Value 




Pr Z 


60 


SirelD(Bcode) 0 


4389 


0.3429 1 


28 


0 


1003 




con temp (Bcode) 3 


3608 


0.7375 4 


56 


< 


0001 




Residual 10 


6090 


0.6789 15 


63 


< 


0001 


65 


Type 3 


Tests 


of Fixed Effects 












Num 


Den 










Effect 


DF 


DF F Value 


Pr > F 




70 


Fingp (Bcode) 


11 


368 16.40 


< 


0001 






Stim 


2 


368 4.68 


0 


0098 






lox (Bcode) 


12 


368 2.18 


0 


0124 
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Table 22 (Continued) 



wt 


1 


368 


8.64 


0.0035 










Estimates 










V 




Standard 










Label 


Estimate 


Error 


DF 


t Value 


Pr > |t| 


Additive Test 11-22 


-0.7500 


0. 8582 


368 


-0.87 


0 


3827 


11-12 


-0.6732 


0.5998 


368 


-1.12 


0. 


2625 


12-22 


-0.07683 


0.7815 


368 


-0.10 


0. 


9217' 


Dominance Test 


0.5963 


1.0975 


368 


0. 54 


0. 


5872 


Recessive Test 


0.8268 


1.5279 


368 


0.54 


0 


5887 


Over Dominance. Test 


-1.4232 


1.2576 


368 


-1.13 


0. 


2585 



Contrasts 

20 

Num t Den 



Label 


DF 


DF 


F Value 


Pr > F 


Additive Test 11-22 




368 


0.76 


0.3827 


11-12 




368 


1.26 


0.2625 


12-22 




368 


0.01 


0.9217 


Dominance Test 




368 


0.30 


0.5872 


Recessive Test 




368 


0.29 


0.5887 


OverDominance Tesc 




368 


1 .28 


0.2585 



30 
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Table 23 



The Mixed Procedure 



Covariance Parameter Estimates 



Cov Parm 



Estimate 



Standard 
Error 



Z 

Value 



10 



15 



20 



25 



30 



35 



SirelD(Bcode) 


0.005452 


0.004008 


1 . 


.36 


0 


contemp{Bcooe) 


0.03915 


0.008774 


4 . 


.46 


< 


Residual 


0.09105 


0.006426 


14 


.17 


< 




Type 3 Tests 


of Fixed 


Effects 








Num 


Den 










Effect 


DF 


DF 


F 


Value 


Pr > F 


Fingp(Bcode) 10 


368 




10.46 


< . 


0001 


Stim 


2 


368 




5.89 


0. 


0030 


lox(Bcode) 


10 


368 




2.28 


0. 


0132 


wt 


1 


368 




7.14 


0. 


0079 




Estimates 














Standard 










Label 


Estimate 


Error 




DF 


t 


Value 


Additive Test 11-22 


0.06263 


0.07965 




368 




0.79 


11-12 


0.09181 


0.05618 




368 




1.63 


12-22 


-0.02918 


0.07277 




368 




-0.40 


Dominance Test 


-0.1210 


0.1028 




368 




-1.18 


Recessive Test 


-0.03345 


0.1419 




368 




-0.24 


OverDominance Test 


0.1544 


0.1171 




368 




1.32 



Pr Z 



.0001 
.0001 



Pr 



|t| 

0.4322 
0. 1031 
0.6887 
0.2398 
0.8137 
0.1879 



40 



Contrasts 
Num Den 

Label DF DF F Value Pr > F 



Additive Test 11-22 1 368 0.62 0.4322 

45 11-12 1 368 2.67 0.1031 

12-22 1 368 0.16 0.6887 

Dominance Test 1 368 1.39 0.2398 

Recessive Test 1 368 0.06 0.8137 

OverDominance Test 1 368 1.74 0.1879 
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Claims ; 

1. A method for assessing the tenderness of meat from an 
animal, comprising the step of testing the animal for the 

5 presence or absence of a genetic marker selected from the 
group consisting of: 

(1) an allele of the gene encoding calpastatin 
(CAST) associated with peak-force variation or 
genetic variation located other than in the CAST 

10 gene which shows allelic association with the 

CAST allele; and 

(2) an allele of the gene encoding lysyl oxidase 
(LOX) associated with variation in instron 
compression of the semitendinosis muscle or 

15 genetic variation located other than in the LOX 

gene which shows allelic association with the 
LOX allele. 

2. A method as claimed in claim 1 wherein the genetic 

2 0 marker is located in the 3' UTR of the CAST gene. 

3. A method as claimed in claim 2 wherein the allele is 
CAST3 D/E allele 1. 

25 4. A method as claimed in claim 2 wherein the allele is 
CAST 3 D/E allele 2 or allele 3. 

5. A method as claimed in claim 1 wherein the genetic 
marker is located in the 5 7 UTR of the CAST gene, 

30 

6. A method as claimed in claim 5 wherein the allele is 
CAST5 allele 3. 

7. A method as claimed in claim 1 wherein the allele is 

3 5 CAST 5 allele 7 or allele 9. 

8. A method as claimed in claim 1 wherein the genetic 
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marker is located in the genomic region embracing the 
coding sequence of the CAST gene. 

9. A method as claimed in claim 1 wherein the genetic 
5 marker is a genetic variation located other than in the 
CAST gene which shows allelic association with either or 
both of the CAST 3 D/E polymorphism and the CAST 5 
polymorphism. 

10 10. A method as claimed in any one of claims 1 to 9, 

further comprising the step of testing for the presence or 

absence of one or more additional genetic markers 

t 

assocxated with peak-force variation. 

15 11. A method as claimed in claim 10 further comprising 

the step of testing for the presence of CAST5 allele 7 or 
allele 9 and/or the absence of CAST 5 allele 3 once the 
presence of CAST D/E allele 2 has been established. 

2 0 12. A method as claimed in claim 1 wherein the genetic 

marker is allele 1 of the LOX polymorphism. 

13. A method as claimed in claim 1 wherein the genetic 
marker is allele 2 of the LOX polymorphism. 

25 

14. A genetic marker for meat tenderness in an animal 
which is a polymorphic form of the CAST gene, being the 
CAST 3 D/E polymorphism or the LOX polymorphism. 

3 0 15. A genetic marker as claimed in claim 14 wherein the 

variable portion of the CAST gene comprises the nucleotide 
sequence set forth in SEQ ID N0:1. 

16. A genetic marker as claimed in claim 14 wherein the 
3 5 variable portion of the CAST gene comprises the nucleotide 
sequence set forth in SEQ ID NO: 2. 
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17 . A genetic marker as claimed in claim 14 wherein the' 
variable portion of the CAST gene comprises the nucleotide 
sequence set fortli in SEQ ID NO: 3. 

5 

18. An isolated DNA molecule comprising the nucleotide 
sequence set forth in SEQ ID NO:l. 

19. An isolated DNA molecule comprising the nucleotide 
10 sequence set forth in SEQ ID NO: 2. 

20. An isolated DNA molecule comprising the nucleotide 
sequence set forth in SEQ ID NO: 3. 

15 21. An isolated DNA molecule consisting of the nucleotide 
sequence set forth in SEQ ID NO:l. 

22. An isolated DNA molecule consisting of the nucleotide 
sequence set forth in SEQ ID NO: 2. 

20 

23. An isolated DNA molecule consisting of the nucleotide 
sequence set forth in SEQ ID NO: 3. 

24. A method for selecting an animal likely to yield meat 
25 of improved tenderness, comprising the steps of: 



35 



30 



(2) 



(1) 



testing the animal for the presence of an allele 
of the gene encoding calpastatin (CAST) 
associated with low peak-force or genetic 
variation located other than in the CAST gene 
which shows allelic association with the CAST 
allele, and/or for the presence of an allele of 
the LOX gene associated with low instron 
compression of the semitendinosis muscle or 
genetic variation located other than in the LOX 
gene which shows allelic association with the 
LOX allele; and 

selecting animals which have the CAST and/or LOX 
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allele and/or genetic variation in allelic 
association therewith. 

25. A method as claimed in claim 24 wherein the allele 
5 tested for is CAST 3 D/E allele 2 and animals which are 

homozygous for this allele are selected. 

26. A method as claimed in claim 25 further comprising 
the step of testing for the presence or absence of one or 

10 more additional alleles of the gene encoding calpastatin 
associated with low peak-force. 

27. A method as claimed in claim 26 wherein the presence 
of CAST 5 allele 7 or allele 9 is tested for, and animals 

15 having either of these alleles in addition to CAST 3 D/E 
allele 2 are selected. 

28. A method as claimed in claim 26 wherein the presence 
of CASTS allele 3 is tested for, and animals having this 

20 allele are rejected despite the presence of CAST3 D/E 
allele 2. 

29. A method as claimed in any one of claims 24 to 27, 
further comprising using the selected animal for breeding. 

25 

30. An oligonucleotide probe for amplification of a 
genetic marker associated with peak-force variation, said 
genetic marker being either an allele of the gene encoding 
calpastatin (CAST) or genetic variation located other than 

30 in the CAST gene which shows allelic association with said 
allele. 

31. An oligonucleotide probe as claimed in claim 30 
selected from the group consisting of: 

3 5 castd 5' cat ttg gaa aac gat gcc tea c 3' 

caste 5' tct acg att age age tea aga gga g 3' 
CAST5U1 5 ' - GTAAAGC C GC AC AAAACAC AC CC AGG - 3 ' 
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CAST5D1 5 ' -GTTTCTGGACCCTCTGGATGAGGAAGCGG-3 ' 

32. An oligonucleotide probe for amplification of a 
genetic marker associated with variation in instron 

5 compression of the semitendinosis muscle, the genetic 

marker being either an allele of the gene encoding lysyl 
oxidase (LOX) or genetic variation located other than in 
the LOX gene which shows allelic association with said 
allele . 

10 

33. An oligonucleotide probe as claimed in claim 31 
selected from the group consisting of: 

LOX K5 : 5' tat cac tga tgt caa acc tg 3' 
LOX K6: 5' act cag gca cca aat age tg 3' 

15 

34. A kit for use in assessing the tenderness of meat 
from an animal and/or selecting an animal likely to yield 
meat of improved tenderness, comprising oligonucleotide 
probes for amplification of at least one genetic marker 

2 0 associated with meat tenderness, said genetic marker being 
either an allele of the gene encoding calpastatin (CAST) 
or genetic variation located other than in the CAST gene 
which shows allelic association with said allele, or an 
allele of the LOX gene associated with low instron 

2 5 compression of the semitendinosis muscle or genetic 

variation located other than in the LOX gene which shows 
allelic association with the LOX allele, and means for 
amplifying DNA. 

30 35. A kit as claimed in claim 34 wherein the 

oligonucletide probes are selected from the group 
consisting of: 

castd 5' cat ttg gaa aac gat gec tea c 3' 

caste 5' tct acg att age age tea aga gga g 3 

3 5 CAST5U1 5 ' - GTAAAGC CGC AC AAAAC AC AC C C AGG - 3 ' 

CAST5D1 5 ' - GTTTC TGGAC CC TC TGGATGAGGAAGC GG - 3 ' 
LOX K5 : 5 ' tat cac tga tgt caa acc tg 3 ' 
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LOX K6: 5' act cag gca cca aat age tg 3' 

36. An animal when selected by the method as defined in 
any one of claims ,24 to 27. 

37. The progeny of an animal as defined in claim 36. 

38. Meat from an animal as defined in claim 36. 

10 39. Meat from the progeny of an animal as defined in 
claim 36. 

40. The use of an animal as defined in claim 36 in 
breeding . 

15 

41. The use of the progeny of an animal as defined in 
claim 3 6 in breeding. 
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